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We revisit the approach proposed by Mukamel and coworkers to describe in- 
-*-* teracting excitons through infinite series of composite-boson operators for both, 

H the system Hamiltonian and the exciton commutator — which, in this approach, is 

'^ properly kept different from its elementary boson value. Instead of free electron-hole 

O operators, as used by Mukamel's group, we here work with composite-exciton opera- 

tors which are physically relevant operators for excited semiconductors. This allows 

<^ us to get all terms of these infinite series explicitly, the first terms of each series 

■^ agreeing with the ones obtained by Mukamel's group when written with electron- 

^^ hole pairs. All these terms nicely read in terms of Pauli and interaction scatterings 

^^ of the composite-exciton many-body theory we have recently proposed. However, 

^!^ even if knowledge of these infinite series now allows to tackle A^-body problems, not 

^ just 2-body problems like third order nonlinear susceptibility x j the necessary 

^ handling of these two infinite series makes this approach far more complicated than 

H 

^ the one we have developed and which barely relies on just four commutators. 

PACS number: 71.35.-y 



1 Introduction 

Most particles known as bosons, are composite particles made of even number of fermions. 
Proper treatment of the underlying Pauli exclusion principle between fermionic compo- 
nents of these particles has been a longstanding problem for decades [1] . Because many- 
body theories for quantum particles were, up to our work, valid for elementary particles 
only [2], sophisticated "bosonization" procedures [3,4] have been proposed to replace 
composite bosons by elementary bosons. These elementary bosons then interact through 
effective scatterings constructed on interactions which exist between their fermionic com- 
ponents, but dressed by "appropriate" fermion exchanges [5]. Although quite popular due 
to the fact that they allow calculations on problems otherwise unsolvable through known 
procedures, such bosonizations have an intrinsic major failure linked to the fact that, by 
replacing two free fermions by one boson, we strongly reduce degrees of freedom of the 
system. This shows up through the fact that, while closure relation for N elementary 
bosons is 



1 



I=^,j:Bi---Bl\v){v\B,,...B,,, (1.1) 



with [Bi, Bj] = 6ij, the one for A^ composite bosons made of two free fermions reads [6] 

I = ^2T.Bl--Bl\v){v\B,,...B,,. (1.2) 

The huge prefactor change from A^! to (A^!)^ makes all sum rules for elementary and com- 
posite bosons, based on this closure relation, irretrievably different whatever are effective 
scatterings produced by bosonization procedures. And indeed, through this prefactor 
difference in closure relations, we have explained [6] the factor 1/2 difference in the link 
between lifetime and sum of transition rates that we had found [7] for composite and 
bosonized excitons. 

Besides bosonization, very few other approaches to interacting composite bosons have 
been proposed. In the late 60's, M. Girardeau [8] suggested to introduce a set of "ideal 
atom operators" in addition to fermionic operators for electrons and protons. These 
operators, which are bosonic by construction, represent all bound states of one atom, but 
not its extended states. They are forced into the problem through a so-called Fock-Tani 
unitary transformation which, in an exact way, transforms one exact atom bound state 
into one ideal-atom state. Unfortunately, this nicely simple result does not hold for A^- 



atom states with N > 2, the procedure turning quite compHcated very fast. This is why, 
although not advocated by Girardeau, we can be tempted by using his procedure as a 
bosonization procedure, i.e., by only keeping ideal-atom operators in transformed states 
and transformed Hamiltonian. We have however shown [9] that, with such a reduction, 
the obtained results for a few relevant physical quantities are at odd from the correct ones, 
even for the sign. The idea to add to fermionic operators for electrons and protons, a set 
of bosonic operators for atomic bound states, is in fact rather awkward because fermionic 
operators form a complete set in themselves; so that Girardeau artificially introduces an 
overcomplete set of operators in a problem already complex, this overcompleteness being 
obviously difficult to handle properly. Precise comparison of Girardeau's procedure with 
the composite-boson many-body theory we have constructed, can be found in reference 
[9]. 

Another approach, still currently used [10-12], has been proposed by Mukamel and 
coworkers in the 90 's. It is based on the fully correct idea that the system Hamiltonian, 
when acting on fermion pairs, can be replaced by an infinite series of pair operators. In this 
approach, the fact that pairs of fermions differ from elementary bosons is kept exactly 
through commutators of pair operators which are also written as infinite series. The 
pair-operators used by Mukamel and coworkers are products of free fermion operators. 
However, as these are not physically relevant pair operators for problems dealing with 
excitons, their calculations turn out very complicated. This is probably why they have 
only derived the first term of the Hamiltonian and pair-commutator series. This thus 
makes their results of possible use for problems restricted to two excitons only. And 
indeed, using them, they have successfully calculated [12] the third order susceptibility 
X^^^ which results from interactions of two unabsorbed photons through their coupling to 
two virtual excitons. 

In this paper, we follow Mukamel and coworkers' idea, but with exciton operators BJ 
instead of products of free-electron and free-hole operators aj^ 6]^ , these exciton operators 
being the ones which create one-electron-hole-pair eigenstates of the system Hamiltonian, 

{H-E,)BJ\v) = 0. (1.3) 

Thanks to the closure relation for A^ composite excitons, Eq.(1.2), it is easy to derive all 
terms of the series for the composite-boson commutator and for the system Hamiltonian 



in an exact way. As expected, prefactors in these infinite series read in terms of tlie two 
key parameters of the composite-boson many-body physics, namely Pauh scatterings for 
fermion exchanges in the absence of fermion interaction, and interaction scatterings for 
fermion interactions in the absence of fermion exchange. 

However, even with these two infinite series at hand exphcitly, so that problems dealing 
with many-body effects between A^ excitons could now be tackled, this approach turns 
out to be definitely far more complicated than the composite-boson many-body theory we 
have recently constructed [13,14]. Indeed, in this new theory, calculations dealing with 
many-body effects between any number A^ of excitons simply reduce to performing a set 
of commutations between exciton operators, according to two commutators for fermion 
exchanges (see Eq.(5) in ref. [15] or Eq. (14) in ref. [13]), namely, 

[5„, Bf] = NBf-\5^, - D„,) - N{N - l)Bf-' ^ A (^ ^,) B], , (1.4) 

n 

[I^™,5r] = Ar5f-^5:{A(::.i) +A(-/)}i?t , (1.5) 

n 

and two commutators for fermion interactions (see Eq.(5) in ref. [16 ] or Eq.(13) in 
ref. [13]), namely, 

[i/, Bf] = NBr\EM + V^) + ^^^(^4-- ^ e (™ B'rrM , (1.6) 

m,n 

[V;t, 5f ] = AT^t^-i ^ ^ ^^ f^ ^t^^t . (1.7) 

m,n 

In these equations, D^i is the exciton "deviation-from-boson operator" defined through 
Eq.(1.4) taken for A^ = 1, namely, 

Dmi = Sm,i — [Bm,Bi] , (1.8) 

while Pauli scattering A [rn'^i) of two "in" excitons (i, j) towards two "out" excitons {m,n) 
follows from Eq.(1.5) taken for A^ = 1 (also see Eq.(4) in ref. [17]), namely, 

[Drm, B]] = E {A fr i) + A (" /) } Bi . (1.9) 

n 

In the same way, "creation potential" V/ of exciton i and interaction scattering ^ (m i) 
follow from Eqs.(1.6) and (1.7) taken for A^ = 1 (also see Eq.(3) in ref.[17]), namely, 

[H,BJ] = E,BJ + V,^ , (1.10) 



[V^\B}] = Y.^[y^)BlBi. (1.11) 

m,n 

Let us note that Eq.(1.8) which basically says that particles are not elementary but 
composite bosons, was known for quite a long time [1,18]. Equation (1.10) is more recent. 
It was introduced by one of us in her theory of exciton optical Stark effect [19,20]. On 
the contrary, Eqs.(1.9) and (1.11) which allow to reach the two elementary scatterings of 
two excitons, namely, A f^i) ^^^ ^ [m])^ ^^^ fundamentally new. They are the keys of 
our composite-boson many-body theory [17]. 

In this paper, we are going to use these four commutators to write the system Hamilto- 
nian H and the deviation-from-boson operator Dmi as infinite series of exciton operators 
BJ. This will allow us to generate physically relevant prefactors for these series. They 
are found to read in terms of exciton energies Ei, interaction scattering of two excitons 
^ (mi)) ^^'^ ttis following sum of Pauh scatterings, 

Ami{n,j) = A (^ i) + A (" f) 

^mi{n,j), shown in Fig.l, corresponds to processes in which excitons {i,j) exchange either 
a hole or an electron, excitons {m,i) having same electron in Ag = A, while they have 
same hole in Xh- Due to electron- hole symmetry, it is quite reasonable to find these two 
processes on the same footing, in the Ami{n,j) factor. 

2 Deviation-from-boson operator 

Let us start with deviation-from-boson operator Dmi defined in Eq.(1.8). Since both, Dmi 
and the product of exciton operators BJ^Bi conserve number of pairs, we can look for Dmi 
as 

oo 

Dm^ = T.Di:L (2.1) 

n=l 

where the most general form for D^l acting in the subspace made of states having p > n 
pairs, can be taken as 

-^mi = z^^mi (^1, . . . ,r„; ri, . . . ,rn)B, . . . B,,Br^ . . . Br„ . (2-2) 

1 "' 

{r} 



We get this series by enforcing it to be such that, when acting on any A^-exciton state 
\^Pn) hnear combination of BJ_^ . . . BJ^\v), it gives the same result as the original operator 
Dmi, namely, 

N 

/^™#Jv) = E^™l^^) ' (2-3) 

n=l 

for any A^-pair state \^n)- We are going to derive the various operators D^^l by iteration, 
starting from n = 1, as we now show. 

Before going further, it is of importance to note that, due to carrier exchanges between 
two excitons, we do have (see Eq.(5) in ref. [17]) 

BlB] = -^X{l'^BlBi. (2.4) 

m,n 

This equation, which comes from the two ways to construct two excitons out of two 
electron-hole pairs, shows that A^-exciton states \iPn) for A^ > 2, as well as operators like 
Djl^l for n > 2 can be written in many different ways, these various forms being related 
through the replacement of any B'^B'' by a sum of B'^B'' according to Eq.(2.4): Just as 
BI_^ . . . B]^ \v) states form an overcomplete set for A^-pair states, B^s form an overcomplete 
set of operators. This, in particular, allows us to guess that, among the various possible 
forms of -D^j , ^^e one which has a physically relevant meaning, is most probably the 
simplest one. We will come back to this fundamental indetermination, linked to exciton 
composite nature, at the end of this section. 

2.1 Derivation of D\J^ 

Let us first consider a one-exciton state \ipi)- By inserting closure relation for one-exciton 
subspace , i.e., Eq.(1.2) taken for A^ = 1, in front of this state, we find 

Drm\'4^l) =Y.D^,Bl^\v){v\Br,\llJl) . (2.5) 

ri 

As DmiBl_^\v) = [Dmi, Bl^]\v) since Dmi\v) = which follows from Eq.(1.8) acting on 
vacuum, we get from Eqs. (1.9,12) 



D„^^\^Pl) = Y. ^Ur[,n)Bl,\v){v\BrA^,) , (2.6) 

r[,ri 

where Ami{n,j) is the combination of Pauli scatterings introduced in Eq.(1.12). 



We then note that projector \v){v\ can be removed from this equation since state 
Br-i^\i/ji) has zero pair while identity operator reduces to \v){v\ for such a state. Conse- 
quently, Eq.(2.6) also reads 



Drml^i) = Y. A™(r;,ri)5;,5,J^i) . 



(2.7) 



r ,ri 



Since this equation is valid for any state \^i), we readily find that operator D^^l such that 
Dmili^i) = -Dmj-|'?/'i) can be taken as 



Dl^= ^A„,(r;,ri)fi,U,, , 



(2i 



with Ami{r[,ri) given in Eq.(1.12). This result is the same as the one given by Mukamel 
and coworkers (see Eqs. (11-13) of Ref.[12]). 



2.2 Derivation of D 



(2) 
mi 



We now consider two-exciton state \ip2)- By inserting closure relation, Eq.(1.2), for two- 
exciton subspace, in front of \ip2), we get 



Dmi\'ip2) 



2! 



Y,D^,BlBl\v){v\Br,Br,\^2) . 



(2.9) 



ri,r2 



To go further, we note that, due to Eqs. (1.9) and (1.12), 

DrruBlXM = {[D^i, Bl^] + B^D^i) Bljv) 

= Y.BI' (A™(r',ri)54 + A™(r',r2)5tJ \v) . (2.10) 

r' 

We then insert this result into Eq.(2.9) and relabel bold indices. By noting that projector 
\v){v\ can again be removed from this equation since B.f^B,,.^\%l)2) also has zero pair, we 
end with 



Dmi\'4^2) 



(20^ . 



J2 Armir[,n)Bl,BlBr.,Br,\^2) ■ 



(2.11) 



,ri,r2 



If we now turn to Dj^l acting on \^2), we note that Brj^\4'2) has one pair so that, if we 
insert closure relation for one-pair subspace in front of this state, we find 



i(i)i 



E 

r[,ri 



Di:,>\^2) = T. ^rrn{r[,n)Bl, 



j:Bl\v){^\Br 



L r2 



BrAi^2) . 



(2.12) 



We can again remove projector \v){v\ from this equation since i?r2-Bri|V^2) has zero pair. 
This readily shows that operator dI^I such that Dmi\^2) = [Dn-j + -0^^) \^p2) can be 
taken as 

Dlnl = -l E AUr'i,ri)Bl,BlBr,Br, . (2.13) 

2.3 Derivation of DI^I 

To better grasp how Dj^l can be constructed by iteration, let us calculate one more D^l 
explicitly. We consider three-pair state {ip^) and inject in front of it, closure relation for 
three-pair subspace. This leads to 

Drm\^3) = j:^ E D.,^,BlXAM(^\Br,B.,,Br,\^,) . (2.14) 

To go further, we do like for Eq.(2.10) and use commutator [Dmi,Bl] given in Eq.(1.9). 
This leads to 

D^^BlXAM = E {AUr',n)BlBl^ + A^,ir',r2)Bl^Bl^ 

r' 

+AUr',r,)Bl^Bl}\v) . (2.15) 

We then inject this result into Eq.(2.14), relabel bold indices and remove projector |f)(f |. 
This leads to 

Drmm = 7t|t^ E Armir'l,n)Bl,BlBl^Br.,Br,Br,\^P3) ■ (2.16) 

y'^-) r[,{r} 

We now turn to Dl^ll^p^). Since 5^ iV^s) has two pairs, we get, by using closure relation 
for two-pair subspace, 

.2!. 



^Sl^s) = (^) E ^rn^ir[,n)Bl,^Bl^BlJ{v\Br.,Br,Br.M3) ■ (2.17) 

r'i,{r} 



(2) 

We do the same for Dl^Hip^) in which Brj^Br^l'ips) has one pair. By collecting all terms, 
we see that operator D^^l such that -Dmi iV^s) = (-0^^ -|- -D^i + Dmi) iV^s) can be taken as 



^S = ^ E AUr[,ri)Bl,BlBlBr,Br,Br, . (2.18) 



^ rlM 



2.4 Derivation of D)^! 

The above results lead us to think that operator D^l can be written as 

Dlnl = 7n E ^rm{r[, n)Bl,Bl . . . B^Br^ . . . 5,,5., , (2.19) 

where 7„ is a numerical prefactor which, in spite of its values for n = (1,2,3), does not 
reduce to (— l)"~^/n. 

To determine 7„, we look for the recursion relation it obeys. To get this recursion 
relation, we follow the procedure we have used for n < 3, namely, we insert closure 
relation for iV-pair subspace in front of state \iPn)- This leads to 

^™|^^) = (^) J2DrrnBl^...Bljv){v\Br^...BrAi^N) . (2.20) 

We then calculate DmiBl_^ . . . Bl^\v) using commutator (1.9); we relabel bold indices and 
remove projector |f)(f |. This gives 

Drm\i'N) = 7^ E Km{r[.r,)B\Bl...Bl^B.,,...Br,B.,,\i^^) . (2.21) 

We then turn to D]^l acting on \iI)n) for n < N . Since state 5^ • • • i?r„ I'^Af) has (A^ — n) 
pairs, closure relation for this subspace leads to 

D'nl\^N) = (j^j^\ In E ^Ur[,n)Bl,Bl...Br„Bl^^...Bl^ 

xBr,...Br^^,Br„...BrMr,) . (2.22) 



(n) 
mi 



By inserting these results into Eq.(2.3), it is easy to show that the form Eq.(2.19) for D 
is indeed valid provided that 7„'s are linked by 

-r- = (^ - E ^^J^z^, . (2.23) 

with 7i = 1. From this equation, it is easy to show that the first 7„'s are 

72 = -1/2 

73 = 1/3 

74 = -11/48 

75 = 11/120, (2.24) 

and so on. . . , with 77V going to zero with increasing N. 



2.5 Other forms of D 



(n) 



As said at the beginning of this section, composite-boson operators BJ form an overcom- 
plete set to describe electron-hole pairs. This is why any given operator acting in A^-pair 
subspace with N > 2, when written in terms of these -Bj's, can appear through different 
expressions. Indeed, due to Eq.(2.4), it is possible to rewrite -BJ/-BL in Eq.(2.19) as 



b:,bi,=-e^(:u)b',.bi. 



r-y. 



?t Dt 



(2.25) 



since Bt^Bl = BlBt,. We then note that 



IZ-^y rl) ^ (m 



ri 



( r' 

^2 


\ 

r2 


r'l 


ri 


m 


' ) 



(2.26) 



where, according to Fig.2(a), A3 is just the exchange scattering between three excitons 
(i,ri,r2). This allows us to replace the first factor of D^l in Eq.(2.19) by 



Y.^rm{r[,n)Bl,Bl = - E A™ (:j :0 Bl,,B\ 



(2.27) 



1 ''2 



While prefactor Am,i(^i,Ti) corresponds to carrier exchanges between two excitons (i,ri) 
leading to (m, r'^) with excitons m and i having either same electron or same hole, pref- 
actor Ami W'^r) corresponds to carrier exchanges between excitons (i,ri,r2) leading to 
(m,r",r2), with excitons m and i also having either same electron or same hole (see 
Fig.2(b)). 

If we keep doing this procedure for bI, BI,^ with B', relabelled as Bl, , and so on ... , 

'2 -^ '2 '1 

we end with Dl^l written in a quite symmetrical form, although far more complicated 



than Eq.(2.19), namely. 



/./ 



D 



(n) 



' ^) in / , ^^mi 

{r'},{r} 



B I . . . B^f Br . . . -Dn 

r\ r' 1^ n '1 



(2.28) 



where the prefactor corresponds to carrier exchanges between (ra+1) excitons (i, ri, . . . , r„) 
leading to {m,r[, . . . ,r'^) in which excitons m and i either have same electron or same 
hole (see Fig. 2(c)). 
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3 System Hamiltonian 



Let us now turn to the system Hamiltonian originally written in terms of fermionic opera- 
tors for free electrons and free holes. It contains kinetic electron and hole contributions. It 
also contains Coulomb interaction between electrons, between holes and between electrons 
and holes. It is actually quite easy to write the electron-hole part of this Hamiltonian in 
terms of excitons. Indeed, by using the link between exciton operators and free-electron 
and free-hole operators, namely. 



B\= E 4A(k/.>keK), (3.1) 

4A=E^l(^|ke,k,), (3.2) 

i 

where (k/i,ke|?) is i exciton wave function in momentum space, we readily find electron- 
hole Coulomb interaction as 

Veh = - Yl ^q'^k.+q^kft-q^khOke 
q.ke.kh 

= -Y.BlB, E K,(^|ke + q,k,-q)(k,,ke|j). (3.3) 

i,j q,ke,kh 

On the contrary, this cannot be done for other parts of the Hamiltonian, namely, kinetic 
energy terms in o)a and Wh and electron-electron and hole-hole Coulomb terms in a^a^aa 
and b'^b^bb. Nevertheless, since both operator H and product of exciton operators Bl^Bi, 
conserve the number of electron-hole pairs, it is a priori possible to write H as 

oo 

H=Y^ /J(") , (3.4) 

n=l 

Hin) ^ y: h^''\T[, . . . , r^; n, . . . , T^)B\ . . . B^^Br^ ...Br,, (3.5) 

{r},{r'} 

SO that H^"^^ acts on states having p > n pairs. This series is determined by enforcing 

N 

H\i,^) = Y.H^-^\i,^), (3.6) 

n=l 

for any state \iI!n) having A^ electron-hole pairs. Here again, ij(") for n > 2 is expected to 
have various forms since any B'^B'^ can be replaced by sum of B'^ B\ according to Eq.(2.4). 
To get the various terms of i/(") expansion, we are again going to extensively use closure 
relation (1.2) for A^-pair states. This will allow us to get one of these possible forms of H 
quite easily. 

11 



3.1 Derivation of H^^^ 

To get H'^^\ we insert closure relation for one-pair states in front of \il)i) in H\iIji). This 
leads to 

H\^,) = Y.HBl\v){v\BrM,) . (3.7) 

We first replace HBl^\v) by Er-^Bl^\v) for exciton operators create one-pair eigenstates of 
the system. We then note that Br^\ipi) has zero pair, so that we can remove projector 
\v){v\ from this equation. This leads to 

H\i,,) = Y.Er,BlBrM,) . (3.8) 

n 

Since H^^'>\iIji) must be equal to H\il)i) for any one-pair state j-i/'i), we readily see that 
H^^'> can be identified with 

H^'^ = Y,E,.,Bl_^Br, . (3.9) 

ri 

3.2 Derivation of //(2) 

We now turn to two-pair subspace. By inserting closure relation for two-pair states in 
front of \ip2), we find 

Hm = (^) E HBlXM{v\Br,Br.A^,) . (3.10) 

We then use Eqs. (1.10,11) to find 

HB^Bllv) = {BlH + Er,Bl + Vr[)Bl\v) 

= {Er, + E^MAM + E e cj ::) BiA^.,jv) . (3.11) 

r[,r!2 

If we insert this result into Eq.(3.10), relabel bold indices and remove projector |f )(f |, we 
end with 

^1^2) = I ^ E Er,Bl^BlBr,B,, + \ E e C? n) BI,B^,B,.,bA I^^) • (3.12) 

Let US now turn to H^^^ acting on \^i). We first note that -Bj,ilV^2) has one pair so 
that closure relation for one-pair subspace leads to 

i7(^)|^i) = E Er,Bl^Bl\v){v\Br2BrMi) , (3-13) 



1*1 ,''2 
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in which we can remove projector \v){v\ since -Bra-BnlV^i) has zero pair. 

This readily shows that H^'^\ such that Hip2) = [H^^^ + H^'^n \4'2), can be identified 
with 



ij(2) = -- E Er,Bl^BlBr,Br, + ] E e (:| ::) Bl,Bl,Br,B^, . (3.14) 



4 , f 



3.3 Derivation of H^^^ 

To grasp how series H is constructed, let us calculate one more H^"'^ explicitly. From 
closure relation for 3-pair states, we find 

1x2 



3! 



^1^3) = U, T.HBt.BlBlJv){v\Br,Br,BrA^s) ■ 



(3.15) 



w 



We then use Eqs. (1.10,11) to find 

HBlBiBlM = (Er, + Er, + Er.MA.Bi» +T.BIb\ \i (* ^A B^ 



^ V* ''2/ ''1 ^ V* '"»/ ''2 



(3.16) 



So that, if we insert this result into Eq.(3.15), relabel bold indices and remove projector 
|w)(i;|, we end with 



H\^z) 



Y^ Er^ BI^ BI^ BI^ Br^ Br,, B, 

{r} 



(3!)2 

CI 






ri 



E ^Cr{:OBl,Bl',Bl,Br.,Br,Br 



m 



(3.17) 



where Cf^ = N(N — l)/2 is the number of ways we can choose 2 excitons among A^. This 
makes C| = 3. 

We now turn to (H'^^^ + H^'^^j I^Jj^) that we calculate by injecting closure relations 
for 2-pair states in front of -B^.J-i/'s) and for one-pair states in front of -Bn-B^alV's)- By 
collecting all these results, we find that H'^^^ such that Hl^ps) = (^H^^^ + H^^^ + i/^^)) iV^g) 
can be identified with 

i7(3) = 1 ^ ErM^BlBl.Br,B.,,Br, - ^ E ^ CI n) Bl,Bl,BlB,,Br,B.,, . (3.18) 

{r} r\,r',,,{r} 
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3.4 Derivation of //(^) 

The above results lead us to think that operator i/*^") can be written as 

H^ = ttn X! ^riBl^ ■ ■ ■ Bl^Br„ ■ ■ ■ Bj.^ 
{r} 

+(^n E 0,7)Bl,Bl,Bl...BlB,.^...B,,, (3.19) 

with an = — 2/3„ = 7„ for n > 1, with 7„ being the prefactor appearing in Dmi series (see 
Eq.(2.23)), while (ai = l,/9i = 0) for n = 1. 

In order to show this nicely simple result, we are going to determine the recursion 
relations which exist between q;„'s and between f3nS. For that, we follow the procedure 
we have previously used, namely, we first insert closure relation for the A^-pair states in 
front of \iPn)- This leads to 



^1^^) = 7^ E HBl . . . Bl^ \v) {v\Br, ...Br, \^n) ■ (3.20) 

{r} 



We then calculate H acting on A^ excitons through Eq.(l.lO). This leads to 

HBl . . . Bljv) = {BlH + Er.Bl + VX)BI . . . BIJv) . (3.21) 

By using Eq.(l.ll), we find 

V^ B^ ...B'' \v) = (\V\bU + B^ V^)b^ ...B^ \v) 



= E e (:i ::) b\bi,bi . . . bi^ \v) + b^v^^.^ . . . < \v) (3.22) 

We iterate the procedure to end with 

HBl . . . Bl \v) = (Er., + ... + ErjBl^ . . . Bl \v) 

E e (:| ::) Bl,Bl,Bt^ . . . BUv) + permutations I , (3.23) 



r' r' 
'l''2 



the total number of these ^ terms being the number of ways C^ we can choose among A^, 
the two excitons having direct Coulomb process. 

If we now relabel bold indices and remove projector |f)(f |, we end with 



^1^^) = 7^ E ErM, ■ ■ ■ BlBr, ...Br, I^Pn) 

w 



+ 7^ E ^{;ri::)Bl,Bl,Bl...BlBr,...Br,\i'N). (3.24) 
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We now turn to /j(") acting on I^Pn) and assume that its general form is indeed given 
by Eq.(3.19). Since state Br^ . . . Br-^\ip]^) has (A^ — n) pairs, we get, by injecting closure 
relation for (A^ — n)-pair subspace. 



^'"'l«-«7v3^S 



Bl . . . BlBl^^ . . . BlJv){v\Br, . . . BrM^) . (3.25) 



We then remove projector \v){v\ as usual. By collecting all these results, we find that 

operator iJ*^") defined through Eq.(3.6) has the form (3.19) provided that q;„'s and /3„'s 

are linked by 

N 'i-i a^ , , 

q2 N-1 n 

with ai = 1 and I3i = 0, due to Eq.(3.9), while P2 = 1/4, due to Eq.(3.14). By comparing 
Eqs.(2.23) and (3.26), we readily see that aN = 'Jn- In order to determine Pn, we first 
note that the recursion relation for aN also reads 

JV 1 ^^ an 

{my {{N-iy.y hUN-ny.y 

{my hi^N-nW ^ • ^ 

Since C^ = A^(A^ — l)/2, this equation is nothing but the recursion relation for /^at 
provided that we replace a^ by {—2I3n) for any N > 2. Consequently, we end with 

j^ = aN = -2pN for N>2 , (3.29) 

while cti = 7i = 1 and Pi = 0, in agreement with our original guess. 

4 Discussion 

4.1 Summary of the results 

The above results lead us to write deviation-from-boson operator Dmi of two composite 
excitons defined as 

[B^, Bj] = 5„,,i - Drm , (4.1) 

15 



through an infinite series of exciton-operator products, according to 

D^^ = E [^ (™ l) + ^ C^ D] Bl, f 1 + E InP^ Br , (4.2) 



ra=2 



^ = E4---4.-i%-i---%- (4-3) 

{j} 

A f ^ [) is the Pauh scattering for carrier exchanges between "in" excitons (i, r) leading 
to "out" excitons {m,r'), with excitons {m,i) having same electron. Electron-hole sym- 
metry is restored through the fact that, in the second term of Eq.(4.2), namely, A (™ ^), 
excitons (171,1) have same hole (see Fig.l). 

7„'s are numerical prefactors which obey the recursion relation 

■,, = JL- V -"' (44) 

with 7i = 1; so that 72 = —1/2, 73 = 1/3, and so on. . . , with 7jv going to zero for 
increasing A^. 

In the same way, the system Hamiltonian, when acting on electron-hole-pair states, 
can be written as an infinite series of exciton-operator products, according to 

H = J2ErBl(l + Y.^r.Pn]Br 



n=2 



+1 E e(:;::)5.\5,tQ-E7nP.-iU.A. (4.5) 

Let us again stress that, since there are two ways to form two excitons out of two 
electron-hole pairs, any product B'^B'' can be written as a sum of B'^B'' according to 
Eq.(2.4). Consequently, it is always possible to rewrite sums appearing in Dmi and H in 
various different ways, Eqs. (4.2-5) being the simplest ones. 

4.2 Comparison with Mukamel and coworkers' results 

In their works, Mukamel and coworkers use free-pair operators B^^ = c\j^^d\^^, where c\^^ 
creates electron on site mi while d)^,^ creates hole on site m2- The fact that they use sites 
while we here use momenta (see Eq.(3.2)) is basically unimportant. They however keep 
the possibility for electrons and holes of these free pairs to differ from free Hamiltonian 
eigenstates. This is why they have nondiagonal contributions in the one-body part of 
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their Hamiltonian, 

iJn = V t^^) c' c + Y^ t^^) d^ d (A 6) 



in\,n\ m2,"2 



As these free-pair states are not physically relevant states to describe a set of A^ interacting 
pairs, Mukamel and coworkers only reach the two first terms of H series, namely, H'^^^ and 
H^'^\ and the first term of Dmi series, their results reading already as rather complicated 
(see Eq.(18) of ref. [10]). To make precise link with their work, we are going to recover 
their results from our compact forms. 

As electron- hole states used by Mukamel and coworkers form complete set, we can 
expand exciton operators Bl in terms of electron-hole operators B\, according to 

Bl = Y.Bi{h\r), (4.7) 

n 

where \r) = Bl\v), while \h) = B}^\v) = cl^_^dl^Jv). 

Using this expansion (4.7), we see that the first term of the H series we have obtained, 
also reads 

ff(l) = Y.Er,BlBr, 

ri 

= Y.hn.nBlB^, (4.8) 

m,n 

where prefactor hmn is nothing but 

hmn = J2Bri{m\i"i){ri\n) = {m\H\n) , (4.9) 

ri 

since H\ri) = i^^Jri). If we now introduce the two-body part of the Hamiltonian as 
written in Eq.(16) of ref.[10], namely, 

H - ^ Vl/^^^ rt r^ r- r, + ^ V T/*^^^ d'^ d'^ d ■ d, -^ W ■ , r^ d'^ d ■ r, 

c „ /_^ ^ m-inxj\kimini^3i^ki ^ „ Z_^ '^m2n2i2fc2 "12 "2 J2"'fc2 /_^ '^ ■rnin2nk2^min2^32^ki i 

(4.10) 
we see that, for H = Hq + He with Hq given in Eq.(4.6), prefactor hmn defined in Eq.(4.9) 
splits as 

h = /i(°) - W 

"'mn mn ''^mim2"l"2 

"'mn ''mini ""1-2, "2 ^ ''m2"2 ™-i'"l ' K^'-^^J 
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in agreement with the result obtained by Mukamel and coworkers for the first term of H 
expansion. 

We now turn to H'^^\ In view of Eq.(3.14), iJ^^) gpij^g ^^ ^(2) ^ ^m _^ ^^2)^ v^heie 

H^ depends on exciton energy E.^ while HI depends on exciton scattering S,- Before 

(2) 
going further, let us note that, since H^ contains exciton energy Er, this term, by con- 
struction, contains a part of electron-hole Coulomb interaction, namely, the one acting 

(2) 
inside one exciton. On the other hand, as He reads in terms of direct Coulomb scat- 
tering ^ L? ^^j between the two excitons (ri,r2), it contains Coulomb scattering resulting 
from electron-electron and hole-hole interactions as well as from electron-hole interac- 
tion between excitons ri and r2 (see Fig.3(a)). Since electron-hole interaction already 
appears in the first-order term H^^^ through Wmim2nin2 in hmn (see Eq.(4.11)), these two 

electron-hole contributions of H'^'^^ must somehow cancel, as we now show. 

(2) 
If we symmetrize H^ and write exciton operators in terms of free pairs according to 

Eq.(4.7), we find 
1 



ri,r2 



-| E BlBlB.B, 

m,n,j,k 



^(m|i/|ri)(ri|A;) E(^l^2)(^2|j) 



(4.12) 



Using Eq.(4.9) and orthogonality of free-pair states, we readily find 

He^ = ~i ^ BlBlBjBkihmkSnJ + hnj5m,k) , (4.13) 

m,n,j,k 

with hmn given in Eq.(4.11). 

(2) 

If we now consider the part of Ht coming from Coulomb scattering ^ between excitons, 
we can rewrite it, using again Eq.(4.7), as 

1 

4 



H^' = - y in'"']B\B\Br,B 

4 ^^. . V''i ''1/ ■'! '"2 '^ 



\Y.BlBiBA E e(:?::)(mK)(nK)(r2|j)(n|^). (4.14) 



m,n,j,k ^'11^21^11^2 



The sum over r's is readily obtained from diagrams of Fig.3(b) in terms of interactions 
V^^^ between electrons, V^"^^ between holes and W between electrons and holes. It reduces 
to 

ymimnki^m^M Sn2,J2 - WmmijikiSm^M Sni,n\ + {1^2}. (4-15) 
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If we now collect the two parts of H^'^\ we can rewrite it as 

H^^^= J2 U^njkBlBiB,Bk + Z , (4.16) 

m,n,j,k 

where Umnjk is just the prefactor obtained by Mukamel and coworkers in Eq.(18) of ref.[10]. 
Operator Z contains all electron-hole contributions. Its precise value reads 



Z = Yl Bl^B^BjEk 

m,n,j,k 



. y'^ mim2kik2^n,j ~r '^nin2JiJ2^'m,,k) 



1 

T \'^min2J2ki^m2,k2'^ni,ji + '^nim2fe2jri"mi,fci0n2,J2/ 



(4.17) 



In order for Eq.(4.16) to agree with the expression of H^"^^ obtained by Mukamel and 
coworkers, operator Z must reduce to zero. This is actually true, as shown by noting that 

^m^n^j^k — C„^_^(l^^^C^_^aj^^aj^Cj^ak2Cki — (^mx\'^n2^ni'^m2)\"'J2^ki(^k2)(^ji ■> \^-^°) 

and by exchanging bold indices {in2 <-* n2) and (J2 ^^ k2) in the sums appearing in Z. 
This explicitly shows that electron-hole interaction does not appear in H^"^' as reasonable 
since, due to Eq.(3.3), Veh can be exactly written in terms of BJBj, or BJ^Bn- 

We thus conclude that expressions of H^^^ and iJ*^^-* given by Mukamel and coworkers 
agree with our compact form of H^^\ As the exciton operators we here use are physically 
relevant operators for interacting electron-hole pairs, we have been able to write the 
whole infinite series for if in a compact form, in terms of these operators. Let us however 
stress that, even with this infinite series now known, it is far simpler to calculate H\%l)j~^) 
through the commutators [i/, B] ] and [Vj , -BJ ], given in Eqs. (1.6,7), than through this 
/j(") series, mostly when the state \'4!n) of interest has many identical excitons, as in usual 
physically relevant situations. 

4.3 Possible use of series expansion for H 

The procedure proposed by Mukamel and coworkers is definitely not a bosonization proce- 
dure, since exact deviation-from-boson operators are a priori kept through their expansion 
as a series of pair operators. We can however be tempted by comparing prefactors ob- 
tained in this expansion of the system Hamiltonian H in terms of exciton operators, with 
effective scatterings produced by bosonization. 
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When truncated to its one and two-body terms, the effective Hamiltonian for bosonized 
excitons reads as 

H = Y: E.BJB, + 2 E ^" (™ i) BlBiB^B, , (4. 19) 

i mnij 

with [Bm, BJ] = 5m,i for elementary bosons. We see that the prefactor of the first term 
of H is nothing but the one of H^^\ If we now consider H^"^^ given in Eq.(3.14), we can 
rewrite it as 

H^'^ = 4 E [^ (a n) - (^n + Er,)i6r'^,r. + 5,,,,j] Bl,Bl,Br,Br, . (4.20) 

Due to the two ways to form two excitons out of two electron-hole pairs which lead to 
Eq.(2.4), we get from this equation used for B'^B'' or BB, 

E e(:j::) 4545.5.. = - e rci::)^;,^;,^.,^., 

ri,r2,r'j^,r!2 fi,r2,r[,r2 

= - E rXi::)Bl,Bl,B,2Br,, (4.21) 

ri,r2,rj,r2 

where .^™ and .^°"*, shown in Fig.3(c,d), are defined as 

r Ci ::) = E A (:? ::) e fe -) , (4.22) 

Pl,P2 

r' CI ::) = E e CI ::) a (- -) . (4.23) 



Pl,P2 



This shows that in Eq.(4.20), C, can be replaced by (— ^™) or (— ^°"*), or even by {a^ — 
6f " - c^°"*) with a + b + c=l. 

If we now turn to the E part of H^'^\ the same Eq.(2.4) used for B'^B'^ or BB leads to 

E {Eri + Er2)Bl^Bl^Br2Bri = - E (Eri + -^rJA Q ^' j 5^, 5^^ 5^2 5^, 

= - E (i?.i+i?r^)A CI ::)5i,5;, 5,25,, (4.24) 

so that -E prefactor in H^'^^ gives rise to two-body scattering between excitons in EX. 
This shows that the second term of H expansion can also be written as 

H^'' = l E ^ CI 3 4^4^-^- (4.25) 

ri,r2,r'i,r!2 

3 CI ::) = 2 b^ CI 3 - '^^" (^1 3 - ^^°" CI 3 - (^^^ + ^^^)^ CI 31 ' (4-26) 
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with a + b + c = I. It is however clear that such a S [y ^j cannot be used as an 
effective scattering between two excitons. Indeed, S depends on energy origin through 
exciton energy Ei which includes the band gap in the case of excitons, while it is physically 
irrelevant to have the band gap entering exciton scattering. Even if we drop these spurious 
E\ terms, this S (^.J ^^j has problem since its ^ part differs from the effective scattering 
between bosonized excitons mostly used in the literature, namely ^ (m''i) = 'C(mi) ~ 
,^°"* (mi)) by at least a factor of 1/2, in addition to the fact that effective Hamiltonians 
with such a ^ are not hermitian: Indeed, in order for H to be hermitian, we must have 
e (mi) = ae (mi) - ^r (mi) - c^"* (mi), with a = a* and 6 = c*. With respect to 
hermiticity, let us recall that Eqs. (4.25,26) are written with composite-boson operators, 
not elementary-boson operators Br'. This makes Eq.(4.21) correct, i.e., H^^' hermitian, 
even for a ^ a* and b ^ c*. 

5 Conclusion 

In this paper, we revisit the procedure proposed by Mukamel and coworkers to approach 
interactions between excitons while keeping their composite nature exactly, through in- 
finite series of composite-boson operators for both, the system Hamiltonian and the 
deviation-from-boson operator of these composite bosons. While Mukamel and coworkers 
use free-electron-hole-pair operators, we here use exciton operators which are physically 
relevant operators for problems dealing with excitons. This allows us to write all terms 
of these two infinite series explicitly. They read in terms of exciton energies as well as 
Pauli and interaction scatterings that appear in the composite-boson many-body theory 
we have recently constructed. We show that the first-order terms found by Mukamel and 
coworkers agree with our results. However, the necessary handling of these two infinite 
series for calculations dealing with A^ excitons makes this approach far more complicated 
than the ones based on the many-body theory for composite bosons we have proposed 
and which only relies on four nicely simple commutators. 
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Figure 1: Shiva diagrams for Ami{n,j) defined in Eq.(1.12). A(^-|j, represented by 
diagram (a), is identical to Ag („ i) , represented by diagram (c), in which m and i have the 
same electron. A (™/), represented by diagram (b), is identical to \h (mi)? represented 
by diagram (d), in which m and i have the same hole. 
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Figure 2: (a) Shiva diagram representation of Eq.(2.27), the summation over the bold in- 
dex r' being performed readily, (b) Shiva diagrams for the prefactor Ami [1^ 11) appearing 
in Eq.(2.28). This prefactor corresponds to carrier exchanges between {i,ri,r2) leading 
to {m,r[,r'2), in which the excitons m and i either have the same electron or the same 
hole, (c) Same as (b) for the prefactor appearing in Eq.(2.29). 
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Figure 3: (a) Direct Coulomb scattering ^ \y |^^ j between the "in" excitons (ri, r2) leading 
to the "out" excitons (^i, T2)- This scattering is "direct" in the sense that the electron-hole 
pairs are coupled similarly in the "in" and "out" states, (b) Diagrammatic representation 
of the sum over {r} appearing in Eq.(4.14) and leading to Eq.(4.15). (c) Shiva diagram 



representation for the "in" exchange Coulomb scattering ^" 



r'2 ra 
r' ri 



defined in Eq.(4.22). 



In this scattering, the Coulomb interactions are between the "in" excitons. (d) Shiva 
diagram representation of the "out" exchange Coulomb scattering ,^°"* \y ^^ j defined in 
Eq.(4.23). In this scattering, the Coulomb interactions are between the "out" excitons. 
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